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We present a quantitative analysis of the flavor physics of a SUSY 50(10) model proposed by 
Chang, Masiero, and Murayama, linking b — > s transitions to the large observed atmospheric 
neutrino mixing. We consider B s —B s mixing and r — > fj,*y and discuss their correlation paying 
careful attention to their different dependences on model parameters. Our analysis shows that 
an order-of-magnitude enhancement of the mass difference with respect to the Standard Model 
is possible and can be reconciled with the most stringent bounds on lepton flavor violation. 
The situation for various AB = 1 decays is discussed qualitatively. 



1 Introduction 

Within the Standard Model, all flavor structure (before electroweak symmetry breaking) resides 
in the Yukawa coupling matrices. However, not all of their parameters are physical. For instance, 

C S m^-H q u Y° d Rj + d Rl $d Ri , Y D = U£Y d U D , 

where the hat here and below denotes a diagonal matrix, and the field redefinition dm = U* Di j d' R j 
eliminates Up without affecting any other term in the Lagrangian, leaving it unphysical. 

In the MSSM additional flavor structure appears in the soft-breaking terms, for instance, 

Ajoft D -d R ™ 2 j R d~R, m\ R = Um\JP. 

As a consequence, the quark-squark-gluino vertex d R id R j(j acquires a flavor-changing structure 
(Ud U)ij. To avoid the associated large flavor-changing neutral currents (FCNC), it is common 
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to consider flavor-universal soft-breaking terms, as in minimal supergravity (mSUGRA). While 
the universality in mSUGRA only holds at a high (Planck or GUT) scale, for the right-handed 
down- type squark masses the radiative corrections are small if tan (5 is not too large, as there are 
no contributions due to the large top Yukawa coupling. The same is true of the slepton masses. 

Within SUSY GUTs, quarks and leptons are unified into irreducible symmetry multiplets. 
Radiative corrections due to yt above the GUT scale then affect full GUT multiplets. This can 
generate mass differences and flavor violation in the down-type-squark and the slepton mass 
matrices, and the mixing angles contained in Ud U can be rendered observable^ 

Furthermore, hadronic and leptonic flavor structures are linked in GUTs. The presence of 
large atmospheric and solar mixing angles in the leptonic sector raises the question whether they 
can manifest themselves in the hadronic sector. Chang, Masiero, and Murayama (CMM) have 
proposed an SO (10) model in which the large v^v T mixing angle can affect transitions between 
right-handed b and s quarksP 

2 Renormalization and particle spectrum in the CMM model 

The 5O(10)-symmetric superpotential has the fornJ^l 

W w = ^16^16 10 H + ^-^16 T Y D 16 10' H 45 H + -^-^16 T Y M 16 16 H 16 H . (1) 

Here the vector 16 in generation space comprises the three generations of matter, and Y u is a 
symmetric 3x3 matrix containing the large top Yukawa coupling. The Higgs fields IOh and 
lO^j contain the MSSM Higgses H u and H^, respectively. The dimension-5 terms, suppressed 
by the Planck mass Mpi, involve further Higgs fields 45h and 16h acquiring VEVs V45 and 
vjq of order the scale Miq, where 50(10) is broken to SU(5). Below the GUT scale Mgut> 
there are the usual four Yukawa matrices of the MSSM (including right-handed neutrinos), with 
Y D oc y D jf^ « Y eT and Y u ^ Y v . Note also that generically tan/3 = O(M 10 /M Pl ). The 
last term in Eq. ^ generates small seesaw neutrino masses. Flavor symmetries are then invoked 
which render Y u and Y M simultaneously diagonal. This defines the (C/)-basis. In this basis the 
remaining Yukawa matrix Y D ~ Y E , giving mass to down-type quarks and charged leptons, 
contains all flavor mixing: 



Vd 

Y D = U C * KM y s ) C/pmns. (2) 

V Vh, 

Here Vckm an d t^PMNS are the quark and leptonic mixing matrices and certain diagonal phase 
matrices have been omitted. We drop the subscript on C/pmns below. 

The soft SUSY-breaking terms are assumed universal near the Planck scale. The only large 
Yukawa coupling yt now renormalizes the 50(10) sfermion mass matrix, keeping it diagonal 
in the (C/)-basis but suppressing the masses of the third generation. At the weak scale, the 
right-handed down-type squark and the left-handed slepton mass matrices thus take the form 



m- 



/m 2 r \ /mf 

m* 3 , m 2 ' {U) = mf 

dm I 1 \ 1 

\ m 2 - - A,- / \ mk - A r 

x d R1 d/ \ h I- 



This nonuniversal structure induces flavor-changing couplings once the charged- lepton and down- 
type-quark Yukawa matrices are diagonalized, as discussed in the introduction. From Eq. EJand 
the GUT relation Y D « Y eT ', it follows that the PMNS mixing matrix U governs the coupling 



of right-handed down-type squarks to gluinos and quarks as well as that of left-handed sleptons 
to charginos and neutralinos and leptons. In particular, the large atmospheric mixing angle will 
appear in the mixing of bn and sr. Given the apparent lack of Yukawa unification for the first two 
generations, which due to their smallness are more sensitive to higher-dimensional operators^ 
it may not be justified to expect a complete correspondence of leptonic and hadronic mixing 
angles. Imposing (approximate) Yukawa unification only for the third generation is sufficient to 
obtain |£7d23| ~ I ^33 1 ~ I | ~ \U T s\ ~ l/\/2, suggesting large FCNC involving second-to- 
third-generation transitions. On the other hand, | C/d 13 1 ~ \U e ^\ is unknown and consistent with 
zero, so (for nearly degenerate first two generations of sparticles) no large effects occur in 1 — ► 2 
and 1 — > 3 transitions. (For instance, the B^ — B^ mixing phase remains 2(3, which ensures 
consistency with the standard CKM fit.) For this reason, the phenomenological analysis below 
is restricted to 2 — > 3 processes. 

Besides tan (3, we choose as low-energy input parameters the approximately universal first- 
generation squark mass nig, the gluino mass m g , and a down-squark ^-parameter = A^^/y^. 
We use the MSSM and GUT RGEs to relate them to the Planck-scale parameters mo, ao, and 
m\/2- The low-energy quantities Aj and then are calculable functions of our low-energy 
inputs, with the biggest contributions coming from the running above M\q. As a further input 
we have the (supersymmetric) /U-parameter, which we take as a free low-energy input. 

A renormalization-group study of the top Yukawa coupling shows that predictivity of the 
model requires yt to remain below a certain value, which translates to a lower bound on tan (3 ~ 3, 
depending on the top mass and (mildly) on supersymmetric parameters^ The largest effects 
occur at this "critical" value. We take m^ ole = 178 GeV and set tan/3 to its critical value in the 
phenomenological analysis below. 



3 B s — B s mixing 

The dominant new contributions to B s — B s mixing in the CMM model are 0(a 2 s ) corrections 
from one- loop box diagrams with gluinos and squarks. One obtains the mixing amplitude^ 

c 2 m 2 

Ml2 = 3/2*7 A ? IP*- + °R) < B * Pl\ B s )■ (3) 

327T Z Mb 3 

Here (B s \Ol\ B s ) is the hadronic matrix element of the Standard Model effective operator 
Ol = slT^^l SLj^bL and Cl is due to Standard Model W— t exchange. The only new operator 
generated is the parity reflection of Ol, with the Wilson coefficient 



Cr _ A§ 8ir 2 a 2 



a s (rrig) 



X 2 GlM^mj 



a s (m h ) 



6/23 

(4) 



and the same matrix element as Ol- A* = V t *Vtb and IA3I = | ^^,3 1 1 1 ~ \ denote the relevant 
flavor-mixing parameters in the Standard Model and the CMM model, and is a dimen- 
sionless function of the squark and gluino masses. Note the enhancement of Cr due to the 
large atmospheric mixing and the large a s . The neutral i? s -meson mass difference is given by 
AMb s = 2|Mi2|, while the phase of Myi is responsible for mixing-induced CP violation. Fig. ^ 
(left) shows a contour plot of the modulus of the new physics contribution normalized to the 
Standard Model prediction in the (ra g -, ad) plane for m~ g = 250 GeV. The shaded area in the plot 
is excluded by experimental constraints on the sparticle spectrum. We observe that the Stan- 
dard Model prediction of about 17.2 ps _1 can be exceeded by up to one order of magnitude. The 
effect decreases rapidly with increasing gluino mass. As the phase of A3 is undetermined, there 
is potentially large mixing-induced CP violation in decays such as B s — > tjj(j) and B s — > ij)rf'\ 
This is also visible from Fig. ^ (right) > where the accessible region in the complex Myi plane 




Figure 1: Left: New-physics contribution to B 3 — B 3 mixing, normalized to SM value, for nig — 250 GeV. Right: 
Allowed region in the complex M12 plane, see text for explanation. 




Figure 2: Diagrams contributing to r — » /17. Approximate mass eigenstates of three different masses are involved, 
two for the left-handed sleptons and one right-handed charged slepton. 



is shown in light grey (orange) and the region excluded by the experimental lower bound ^ 
AMb s > 14.5 ps _1 is shown in dark grey (blue). The SM prediction, almost real, is also indi- 
cated. 



4 AF = 1 processes: Constraints and predictions 

4-1 Constraint from BR(t — > ^7) 

In the leptonic sector a strong constraint comes from the process r — > fij. The Standard Model 
contribution is negligible, and the supersymmetric contributions^ are due to the diagrams shown 
in Fig. [21 Note that the amplitude involves a chirality flip, either via the equation of motion on 
the external r line, through a Yukawa-like coupling in the vertex, or through an internal LR flip 
due to left-right mixing in the slepton mass matrix (only for the neutralino contribution) . In the 
CMM model, this LR mixing is small and can be treated in the mass-insertion approximation. 
Such an expansion would not be well justified for the LL mixing due to the large mixing angle 
present there. The branching ratio takes the (schematic) form 

BR(t -> fi-y) oc I Cy 1 2 , C7 = . 2 U t3 U* 3 A [m- v m^o,- , y) . 

(4:TT) Z sin t>w V / 

As an important difference to B s — B s mixing, it depends on the //-parameter, which is un- 
correlated with the squark and slepton masses in the CMM model. Numerically we find that 
the chargino-sneutrino diagrams dominate the amplitude for < 1000 GeV. Furthermore, the 
dominant piece exhibits decoupling as 1/fj, for |/i| — ► 00. The effect of this is visible in Fig. |31 
(left), where the dependence of the branching ratio (keeping all terms) on [i is shown for fixed 
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Figure 3: Left: Dependence of BR(r — > /ry) on fj,. The dashed line corresponds to the experimental upper 
bound. For parameter values besides fj,, see text. Right: Same as Fig. (left), but including contours where the 
experimental limit on BR(r — > /17) is reached. The long-dash-dotted line corresponds to /i = —300 GeV, the 

short-dash-dotted one to /i = —600 GeV. c 

parameters nig = 250GeV,mq = 1000 GeV and for ad = 0.75 m^. We indeed verify from the 
figure that the impact of the BaBar bound^ BR(t — > fij) < 6.8 x 10 8 strongly depends on 
fx. For a given /x, the contour where the bound is reached divides the plane of Fig. ^ (left) 
into an allowed and an excluded region. This is illustrated in Fig. El (right). We conclude that 
without specifying n, large B s — B s mixing cannot be excluded by r — > fij: for instance, for 
[i < —600 GeV the impact of the bound almost disappears. 

4-2 Bounds from AB = 1 processes 

While so far there exist only lower bounds on the B s — B s mass difference, one can also try 
to constrain the CMM model through information on b — > s transitions such as B& — > X s j, 
Bd — > X s l + l~ and B^ — > 4>Ks, which have been observedP^We have so far not studied these 
in comparable detail. However, we can make the following remarks: B^ — > X^, like r — > /ry, 
is sensitive to left-right mixing, which is generally not large in the CMM model. Furthermore, 
the chiral flavor structure of the model would primarily generate the mirror image Q' 7 ^ of the 
SM magnetic penguin operator Q7 7 ; then the SM and CMM contributions do not interfere 
due to the different chirality structureP In view of the still substantial theoretical uncertainties 
in this decay, it is possible that information on this mode poses no additional constraints. 
Similar comments apply to Bj — > X s l + l~ . Nevertheless, a complete analysis should include a 
more detailed treatment. The branching ratios for loop- induced hadronic decay modes such as 
B& ~ ¥ 4>Ks suffer from large theoretical uncertainties, precluding any effective constraint. 

4-3 Mixing- decay interference in b — > s penguin decays 

Unlike the decay rates, the time-dependent CP asymmetries in penguin-dominated B^ — > MK$ 
decays with a vector or pseudoscalar meson M in the final state are theoretically rather clean 
and close to the one in B^ — > J/ipKs in the Standard Model. This is due to the fact that, 



C A color version of this plot is included in the slides for this talk, available on the Moriond website. 



the decay amplitudes being nearly real, all are dominated by the single weak phase 2/3 entering 
through the B^ — B^ mixing amplitude. 

While in the latter mode the decay amplitude arises at tree level, in the former it is loop- 
induced and both its modulus and weak phase are sensitive to new physics, destroying the 
relationP^ Indeed there is a well-known discrepancy between the experimental world average 
over penguin modes, —r]cpSf(b — > s penguin) = 0.43 ± 0.07, and that over charmonium modes, 
sin2/?(6 -> ccs) = 0.726 ± 0.037P] 

While we have not yet performed a detailed numerical analysis, analytical results obtained 
within QCD factorization for the relevant operator matrix elements show that the parity re- 
flection Q' 8g of the SM chromomagnetic penguin operator could give an 0(1) correction to the 
Standard Model amplitude, again with an unconstrained phase 6 correlated with the B s — B s 
mixing phase. We remark that in the limit where this contribution dominates the amplitude, 
the CP asymmetry would measure sin(2/3 + 26) in all cases. In the general cas e of interference 
with th e St andard Model contribution, there is dependence on the final stately A recent generic 
analysis^l relevant to the CMM model suggests, however, that right-handed squark mixing can 
account for the data in both vector and pseudoscalar final states. 



5 Conclusion 



We have studied a SUSY SO (10) model implying correlations between neutrino mixing angles 
and the flavor structure of certain sfermion mass matrices. Large effects are found possible 
for B s — B s mixing and r — ► fj,j. Both are correlated, but the additional dependence of the 
latter on the //-parameter allows to satisfy the experimental bound while retaining a large mass 
difference and/or CP violation in the former. A large chromomagnetic penguin in b — > s penguin 
decays is suggested within QCD factorization. Whether this can generate time-dependent CP 
asymmetries at the observed central values remains to be resolved through a numerical analysis. 
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